How orbital angular momentum affects beam shifts in optical reflection 
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It is well known that reflection of a Gaussian light beam (TEMoo) by a planar dielectric interface 
leads to four beam shifts when compared to the geometrical-optics prediction. These are the spatial 
Goos-Hanchen (GH) shift, the angular GH shift, the spatial Imbert-Fedorov (IF) shift and the 
angular IF shift. We report here, theoretically and experimentally, that endowing the beam with 
Orbital Angular Momentum (OAM) leads to coupling of these four shifts; this is described by a 
4x4 mixing matrix. 

PACS numbers: 42.79.-e, 41.20.Jb, 42.25.Gy, 78.20.-e 
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I. INTRODUCTION 

The reflection of a light beam by a mirror shows sub- 
tle aspects that were first conjectured by Newton 
the center of the reflected beam may show a small spa- 
tial shift in the plane of incidence, relative to the po- 
sition predicted by geometrical optics. This shift has 
been named after Goos and Hanchen (GH) who were the 
first to observe it in total internal reflection (TIR) |2j. 
Additionally, there is a spatial shift perpendicular to the 
lane of incidence, the so-called Imbert-Fedorov shift (IF) 
H|. There exist also angular GH and IF shifts, both 
of which have been demonstrated recently in external re- 
flection 0, [g| . The angular shifts can be seen as shifts in 
wave- vector space [6|-[8j. All these shifts depend on the 
polarization of the incident photons. Accurate calcula- 
tions of either GH or IF shifts (or both) can be found 
in Refs. [ol-Hlj. In more recent years the GH shift has 
been studied in a large diversity of cases, ranging from 
photonic crystals [12] to neutron optics [13| . 

We are interested in the question how these beam shifts 
are affected when the light beam is endowed with Orbital 
Angular Momentum (OAM). OAM is a relatively novel 
degree of freedom of a light beam that has found ap- 
plications from optical tweezers to quantum information 
science [H], [H| . Theoretically, a treatment of the effect 
of OAM on beam shifts has aready been given, first by 
Fedosoyev [l6l-[l8j and then by Bliokh et al. [19]. Here we 
prefer to develop our own theoretical treatment based on 
straightforward application of Snell's law and the Fresnel 
equations, in order to derive a unified matrix formalism 
for the four basic shifts: spatial GH, angular GH, spatial 
IF and angular IF. Experimentally, Okuda and Sasada 
have studied the deformation of an OAM carrying beam 
by TIR very close to the critical angle [20]; however, they 
did not report GH and IF shifts. Dasgupta and Gupta 
have measured the IF shift of an OAM beam reflected by 
a dielectric interface, but only for the spatial case [21] . 

It is the purpose of this article to report a theoretical 
and experimental study of the effect of OAM on the four 
basic shifts: spatial GH, angular GH, spatial IF and an- 



gular IF. We find that these shifts are coupled by OAM; 
this is described by an OAM dependent 4x4 mixing 
matrix. We have experimentally confirmed this mixed 
occurrence of GH and IF shifts. 



II. THEORY 

In this section we furnish a thorough theoretical anal- 
ysis for the problem of the reflection of an OAM carrying 
light beam by a dielectric interface. 

Consider a monochromatic beam containing a continu- 
ous distribution of wave- vectors k centered around ko = 
hoZi, where Z{ is a unit vector along the central propa- 
gation direction of the incident beam: k = /cok = ko + q, 
with q = qr + <7L^i and q\t • z% =0. Using the notation 
of Fig. 1, we write qr/ko = sin a and qr/ko — ^ — cos a 
with qr = |qr| and a = arccos(k-ko). A collimated beam 
has a narrow distribution of wave- vectors around ko such 
that sin a = a <C 1 with qr/ko = a <C 1 and qr/ko — 
(<2tAo) 2 / 2 - Thus > if we write k = k (xiU + yV + ZiW) 
with W — y/l — U 2 — V 2 , we can assume U, V <C 1 with- 




FIG. 1: (Color online) Geometry of beam reflection at a di- 
electric interface. The reflected wavevector k is the mirror 
image of the incident wavevector k. Inset shows quadrant 
detector with sensitive areas a, 6, c, d. 
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out significant error. Let E J (r,£) be the electric field of 
the incident beam. Upon reflection this field evolves to 
E jR (r, t) which is to be found. From the linearity of the 
wave equation it follows that E jR (r,t) can be determined 
by studying the action of the interface upon each plane 
wave field 

A\k) = /j_(k) exp (ik • r - iut) (1) 

that constitutes E J (r, £), with uj = |k|c. According to 
Refs. [H, [23|, we assume the polarization dependent 
amplitude of A 1 (k) equal to /j_(k) = f — k(k • /) = 
a p (k) £ k + a s (k) y k with 

a p {k)=f p + Vf s cotO, a s (k) = f s -Vf p cotO, (2) 

up to first order in J7, V, and # = arccos(£^ • z). Here 
/ = + fsV is a unit complex vector that fixes the 
polarization of the incident beam, and y v = z x v/\z x 
v|, x w = y w x v denote a pair of mutually orthogonal 
real unit vectors that together with the arbitrary vector 
v = v/|v| form a right-handed Cartesian reference frame 
K w = {cCv,2/v,v} attached to v. 

When the beam is reflected at the interface, each plane 
wave evolves as: A 1 (k) — » A R (k) where 

A R (k) = [r p (0k)a p x^ + r s (0k)a s y^] x (r,t), (3) 

and x(f, t) = exp (ik • r — iut) = exp (ik • f — iut). The 
notation v indicates the mirror image of the vector v with 
respect to the interface: v = v — 2z(z • v), with v • u = 
v • u [24]. Moreover, r p {0\^) and r s (#k) are the Fresnel 
reflection coefficients at incidence angle #k = arccos(k-i) 
for p and s waves, respectively. By direct calculation it 
is not difficult to show that, up to first order in C7, V, 

x^. = x r — V cot Oy + U z r , (4) 
y^ = V cot 0x r + y — Vz r , (5) 
rx^) = r x + Ur' x , (6) 

where A G {p, s}, r A = r\{6), and = dr\(0)/dO. With 
the use of Eqs. ([211 16]) into Eq. © we obtain 

A*(k) = x r A R (k) + yAf (k) + £ r Af (k), (7) 

where, up to first order in U, V, 

Af(k) = f x r x (1 + iX x U - iY x V) x (f , t), (8) 
4£(k) = (/ p r p f/ - f s r s V) X (r, t). (9) 

Here, we have defined 

x 5 = -^, + (10) 

with X s = XJ and Y s = -YJ ^ . The limit of 
specular reflection is achieved by letting r p — » 1 and 
r s — >> — 1 where Eq. (pTU]) reduces to X p = = X s 



and Yp = = Notice that from Eqs. flEJE]) it fol- 
lows that for a paraxial beam the longitudinal electric 
field energy density \A R | 2 scales as ~ a 2 and it is there- 
fore negligible with respect to the transverse electric field 
energy density \A R \ 2 + \A R \ 2 that scales as ~ 1 + 2a. 
Thus, up to first order in a, we can neglect the longi- 
tudinal term A% (k) and write A R (k) ~ £ r -4p + . 
Moreover, for small shifts X A and Y A one can write 
1 + iX x U - iY x V ~ exp(iX x U - iY x V), and in the 
Cartesian coordinate system attached to the reflected 
beam x(r , t) = ex.p[i(-UX + VY+WZ)} exp(-iut), with 
X = k^x ri Y = koy, and Z = k^z ri where z r is the dis- 
tance from the waist of the incident beam to the quad- 
rant detector measured along the trajectory of the beam. 
Thus, Eq. dSj) can be rewritten as: 

A R {k) ~ fxrxx(-X + X X , Y - Y A , Z, t). (11) 

The passage from the single plane wave field A R (k) 
to the total electric field E jR (r, t) is realized by sub- 
stituting the plane wave scalar amplitude x( r ^) m ^o 
Eq. {UJ, with the electric field scalar amplitude 
E(r, t) describing the spatial distribution of the in- 
cident beam. In the present case, as we want to 
study the behavior under reflection of OAM beams, 
we choose E(r,t) = ip£(r)exp(—iu)t), being ^( r ) the 
Laguerre- Gauss paraxial field with OAM index i G 
{0, ±1, ±2, . . .} and radial index p = 0: i/j£(X, Y, Z) oc 
exp[-(X 2 + Y 2 )/(2A + i2Z)]{X + is l Y) }(l \ with s £ = 
sign(£) and A = ko(kow 2 ] /2) denoting the dimensionless 
Rayleigh range of the beam with waist wo [25[. Thus, 
the transverse electric field of a Laguerre- Gauss beam 
reflected by a plane interface can be written as: 

E?(r, t) ~ fxrxM-X + X X ,Y - Y x , Z) exp(-iwt). 

(12) 

In this expression the terms X\ and Y\ are responsi- 
ble for the GH and IF shifts of the center of 
the beam, respectively. These displacements can be as- 
sessed by measuring the position of the center of the 
reflected beam with a quadrant detector centered at 
x r = 0, y = along the reference axis z r attached to 
the reflected central wave vector ko = koz r . A quadrant 
detector has four sensitive areas, denoted with a, 6, c, d 
in the inset of Fig. 1, each delivering a photocurrent 
7 a , ifo, 7 C , Id respectively, when illuminated. The two cur- 
rents I x = {I a +h)-(Ic+Id) and/^ = (h+Id)-{Ia+I c ) 
are thus proportional to the x- and the ^/-displacement 
of the beam intensity distribution relative to the center 
of the detector, respectively. 

Hi — , ipo ( — X + X\ , Y — Y\ ) reduces to a shifted fun- 
damental Gaussian beam, and in the hypothesis of small 
deviations X\ , Y\ <C 1 , a straightforward calculation fur- 
nishes 

y = N (a gh + |@ghJ , T -f = N (a if + jQwJ , 

(13) 
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where / = I a + 1^ + I c + i^, and A^o = \/2/(ira 2 ) with 
cr 2 = (A/2) ^1 + Z 2 /A 2 . Here we have defined the two 
spatial (A) and the two angular (0) shifts 

A GH = Yl W *MX\), Aif = w ARe(F A ), (14) 

X=p,s X=p,s 

and 




e GH = ™Alm(X A ), 9if = ^Alm(y A ), (15) 

\=p,s \=p,s 

respectively, where the nonnegative coefficients w\ are 
defined as the fraction of the electric field energy with 
polarization A = p, s in the reflected beam: 



\rxfx\ 2 



\r P f P \ 2 + \r s f s \^ 
If t ^ 0, Eq. (HSJ) becomes 



(16) 



TV, ( Af F + -ef F j , 

(17) 
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where N £ = N T(\l\ + 1/2)/ [r(|/| + 1)^ (T(x) denotes 
the Gamma function), and 



• (18) 



Equation ([15]) clearly displays the mixing between spatial 
and angular GH and IF shift, occurring only for i ^ 
0, and it is in agreement with the results presented in 
Ref. [19], apart from the factor "2" in front of £ [26|. 
Notice that the polarization dependence of the four in- 
dependent shifts on the left side of Eq. JT5J) resides in 
the four -^-independent shifts on the right side of the same 
equation. It turns out that the 4x4 mixing matrix itself 
is polarization-independent. It should be noticed that 
in TIR, contrary to partial reflection, both GH and IF 
angular shifts Oqh and Oif are identically zero since the 
Fresnel coefficients are purely imaginary [27j. Thus, in 
this case it follows from Eq. ([T8]) that mixing vanishes. 



III. EXPERIMENTAL SET-UP 

Our experimental set-up is shown in Fig. 2. A home 
built HeNe laser (Ao = 633 nm) is forced to operate in 
a single higher-order Hermite- Gaussian (HG nm ) mode 
with m = by insertion of a 40-/im diameter wire nor- 
mal to the axis of the laser cavity [28]. The HG n o beam 
is sent through an astigmatic mode converter consisting 
of two cylinder lenses, with their common axis oriented 
at 45° relative to the intra-cavity wire. This introduces 
a Gouy phase which converts the HG n o beam in a LG^ P 



FIG. 2: (Color online) Experimental set-up. The insets show 
the HGio and LGio mode profiles. The quadrant detector 
measures the OAM controlled shift of the reflected beam in 
the plane of incidence (GH shift) and perpendicular to it (IF 
shift). Both the GH and the IF shift have a spatial and an 
angular contribution. See text for further details. 



beam with £ = n and p = [28]. Lenses 1 and 2 are 
used for mode matching; the beam leaving lens 2 is col- 
limated with a waist parameter wo = 775 /im, a power 
of typically 600 fiW and a polarization set by a linear 
polarizer. We have incorporated the option to greatly 
enhance the angular spread of the beam by inserting lens 
3 (/ = 70 mm), leading to wo = 19 /im. Either with or 
without lens 3 present, the beam is externally reflected 
by the base plane of a glass prism (BK7, n = 1.51). 
We measured the polarization-differential shifts of the re- 
flected LG^o beam with a calibrated quadrant detector. 
We also obtained these shifts for the fundamental LGoo 
beam (= TEMoo) by simply removing the intra-cavity 
wire from the HeNe laser. 

It follows from Eqs. (|13H17j) that using a collimated in- 
cident beam, i.e. A ^> Z, leads to total predominance of 
the spatial shift. On the other hand, the use of a focused 
beam, i.e. A<Z, leads to total predominance of the an- 
gular shift. These two extreme cases were realized in our 
experiment by removal respectively insertion of lens 3. 
Specifically, the value of the Rayleigh range L — k^w^j2 
was 2.96 m and 1.8 mm, respectively; as standard we 
have chosen the distance z r between the beam waist and 
the quadrant detector to be 9.5 cm. We experimentally 
checked the angular nature of the shift (where expected) 
by verifying that the detector signal depended linearly 
on changes in z r . 

We performed all measurements by periodically 
(2.5 Hz) switching the polarization of the incident 
beam with a liquid-crystal variable retarder and by syn- 
chronously measuring (with a lock-in amplifier) the rel- 
ative beamposition for one polarization with respect to 
the other 0, [29[ . Experimentally we were restricted to 
using the first-order LG modes {I = ±1) by the low gain 
of the HeNe laser. 



IV. EXPERIMENTAL RESULTS AND 
COMPARISON WITH THEORY 

Our experimental results for the polarization- 
differential shifts versus the angle of incidence are re- 
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ported in Fig. 3, together with the theoretical curves 
(£ — and £ = ±1) which are based upon Eqs. ([13JJ17p . 
The four panels show the spatial and angular varieties of 




(degrees) 9 (degrees) 



FIG. 3: (Color online) Reflective beam shift for partial di- 
electric reflection from an air-glass interface as a function 
of the angle of incidence. Plotted curves are the theoretical 
polarization-differential shifts for the two polarizations indi- 
cated in each panel. Experimental data and theoretical curves 
refer to t — and £ = ±1. The panels display the spatial GH 
shift (a), angular IF shift (b), spatial IF shift (c) and angular 
GH shift (d). Here ko — 27r/Ao; see text for further details. 

GH and IF shifts. Note that we have plotted here the 
true GH and IF shifts A/fcrj an d 0/ A, respectively, and 
not the dimensionless shifts A and 0. In each individual 
case the polarization modulation basis has been chosen 
such as to maximize the magnitude of the OAM effect. 
Fig. 4 shows a cartoon-like representation of the four 
cases that we address. The overall agreement between 
experiment and theory is reasonable if we realize that 
there is no fitting parameter involved; we ascribe the re- 
maining discrepancies to insufficient modal purity of the 
LGio beam (we are very sensitive to this since we use a 
quadrant detector). 

Fig. 3a shows the spatial GH shift for a polarization 
basis of diagonal linear polarizations. In this case, the 
GH shift is absent for £ = but it appears for £ = ±1; 
the sign of the shift reverses when going from £ = +1 to 
£ = — 1. In Fig. 3b we show that the angular IF shift is 
different for £ = and £ = ±1, using again diagonal linear 
polarizations. No difference occurs for £ = +1 versus £ = 
— 1. Proceeding to Fig. 3c we observe an angular GH shift 
when using a linear polarization basis for both £ = 

and £ = ±1. Both cases show a dispersive resonance at 
the Brewster angle; for £ = these experimental results 
have been reported recently [Bj] whereas the data for £ = 
±1 (with opposite sign for £ = +1 and £ = —1) are 



(a) (b) 




FIG. 4: (Color online) Reflection of an orbital-angular- 
momentum (OAM) carrying Laguerre- Gaussian (LG) beam 
at a dielectric interface. Depending on the input polarization 
this may lead to a longitudinal shift (Goos-Hanchen effect) 
or to a transverse shift (Imbert-Fedorov effect), where longi- 
tudinal and transverse refer to the plane of incidence. Each 
of these shifts consists of a spatial part and an angular part, 
which are observed, respectively, in the near field and in the 
far field of the LG beam. The magnitude of the shifts in- 
creases with the OAM index £. Panel (a) shows the spatial 
Goos-Hanchen effect, panel (b) the angular Imbert-Fedorov 
effect, panel (c) the spatial Imbert-Fedorov effect and panel 
(d) the angular Goos-Hanchen effect. Note that A and O 
are dimensionless quantities; A and ko are the dimensionless 
Rayleigh range and the wavenumber of the LG beam, respec- 
tively. See text for further details. 



new. Fig. 3d shows the OAM dependence of the spatial 
IF shift, observed in a linear polarization basis (p,p — 2°) 
[2l|, [30| . Here the shift is zero for £ = whereas it shows a 
dispersive Brewster resonance for £ = ±1 (with opposite 
sign for £ = +1 and £ — —1). 

Finally, we have confirmed experimentally that OAM 
did not affect spatial and angular GH and IF shifts in 
the TIR case (not shown); TIR was realized by flipping 
the glass prism in Fig. 2. 



V. CONCLUSIONS 

We have presented a unified theoretical description of 
how the Orbital Angular Momentum (OAM) of a light 
beam affects its kinematic degrees of freedom when the 
beam is reflected by a dielectric interface. Without OAM 
the reflection leads to four beam shifts relative to geomet- 
rical optics, namely the Goos-Hanchen (GH) and Imbert- 
Fedorov (IF) shifts, each of which may have a positional 
and an angular part. We introduce a 4 x 4 polarization- 
independent (but ^-dependent) coupling matrix that de- 
scribes the OAM induced mixing of these four shifts when 
using a quadrant detector. Experimentally, we have con- 
firmed this theory by measuring the four shifts as a func- 
tion of the angle of incidence, for OAM values £ = 
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and ±1. We have observed for the first time the OAM 
induced spatial GH shift as well as the OAM affected an- 
gular GH and IF shifts (see Fig. 3 a, b, c). Extension of 
all this from reflection to transmission (i.e. refraction) is 
straightforward. 

Understanding these effects is important since they 
generally affect control of OAM beams by mirrors and 
lenses. The angular shifts are particularly interesting 
from a metrology point of view, both classically and 
quantum mechanically, since the corresponding trans- 
verse excursion of the beam center grows without lim- 
its when the beam propagates; this greatly promotes its 
detectability [!,[3l|. 
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